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We compute a full set of the first order QED corrections to the parity-violating observables in 
polarized M0ller scattering. We employ a covariant method of removing infrared divergences, com- 
puting corrections without introducing any unphysical parameters. When applied to the kinematics 
of the SLAC E158 experiment, the QED corrections reduce the parity violating asymmetry by 4.5%. 
(^■^ ' We combine our results with the previous calculations of the first-order electroweak corrections and 

\ obtain the complete (0)(a) prescription for relating the experimental asymmetry Alr to the low- 

^vQ , energy value of the weak mixing angle sin'^ 6'vi/. Our results are applicable to the recent measurement 

of Alr by the SLAC E158 collaboration, as well as to the future parity violation experiments. 
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I. INTRODUCTION 



Measurements of the parity- violating observables in electron scattering provide information about low-energy struc- 
ture of weak neutral current processes. Such observables arise from the interference between the weak and electro- 
magnetic amplitudes and are sensitive to the electroweak couplings. In the Standard Model, the couplings of 
QQ . the boson to the fermions are determined by the weak mixing angle 6*-^, which has been measured with high 
^sj ■ precision at the Z'^ resonance Precision measurements of the low-energy parity- violating observables provide an 
independent determination of the weak mixing angle, directly test higher order electroweak corrections, and provide 
' strong independent constraints on the new physics contributions at the TeV scales Q . 
J;0 \ The experiment E158 at the Stanford Linear Accelerator Center (SLAC) has measured the parity-violating left- 
. right asymmetry Atr in M0ller scattering of polarized 50 GeV electrons off unpolarized atomic electrons in a liquid 

■ hydrogen target The final uncertainty on Alr is about 10%. The measurement of Alr translates into the 
O-ij measurement of s\t? 9w with a precision of (7{sit? 9w) ~ 0.001 at low momentum transfer w 0.03 GeV^/c^, and 
O is sensitive to both electroweak one-loop radiative corrections and new physics phenomena at the TeV scales. 

1^ , A precise comparison of the experimental results with the Standard Model predictions requires detailed under- 

■ standing of the radiative corrections, including effects of the QED bremsstrahlung. The leading-order electroweak 
\ corrections to the M0ller scattering have been computed a number of years ago 0, [g . The authors of Ref. factor- 

• , ized out the soft bremsstrahlung contribution, but did not include the effects of hard bremsstrahlung, arguing that 

■ their effects are small, and would require the knowledge of experimental kinematics and acceptance. In this paper, 
improving on the calculation of Ref. Q, we present a complete calculation of the first order QED corrections. The 
detailed analysis shows that the QED corrections are indeed small but not insignificant, compared to the systematic 
uncertainties of E158 and the projected uncertainties of the future parity violation experiments. 

Experimentally, M0ller scattering is often used to measure polarization of electron beams 8]. In such measurements, 
both beam and target electrons are polarized, and electroweak effects can typicall y b e neglected. QED corrections 
to the parity-conserving polarized M0ller scattering have been computed in Ref. lOlllOj and are relatively important, 
compared to the typical precision of M0ller polarimeters. Similar to Ref. Il0j| . we perform our calculations in 
the covariant framework of Bardin and Shumeiko pdj . which allows us to cancel out infrared divergences without 
introducing unphysical parameters (such as a frame-dependent cutoff l^E that separates the soft bremsstrahlung 
region from the hard bremsstrahlung contributions). 

This paper is organized as follows. Section ^1 introduces the kinematics of the M0ller scattering and Born cross 
section and parity-violating asymmetry. Section IlIII explains the regularization of infrared divergences. Section IVII 
presents numerical results applied to the kinematics of SLAC El 58 experiment. 



II. DEFINITIONS, BORN CROSS SECTION, AND KINEMATICS 



The lowest-order (Born) cross section for M0ller scattering, defined by the Feynnian diagrams in Fig. ^ can be 
written as 
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FIG. 1: Born-level diagrams for M0ller scattering: t-channel (left) and it-channel (right). 



W ^ [X'l(u^D\t)D^{t)+t^D\u)D^{u)) + )C_ls^{D'{t) + D\t)){D^t)+D^{u))], (1) 

where the four-momenta of the initial and final electrons fci, pi and fc2, P2 (see Fig^ are combined to form Mandel- 
stam invariants 

S - (fci + pi)2, i - (fci - fc2)^ U - (fc2 - Pi)' . (2) 

In Eq. Q and hereafter, we use a short-hand notation cr = da/dy. 

The Born cross section in Eq. Q is written in terms of the photon and propagators 

D\k)^D^'^^-^ i^^J,Z) (3) 

and the coupling factors 

X±. — Ai^Aiy ± X2^b^2t ■ (4) 
The latter in turn depend on the polarizations of the beam (pb) and target (pr) electrons: 

'^lB(T) ^ -PB{T)X^A ' ^2B(r) = -PB{T)Xv^ (5) 

A'^^ = w'vJ + a'a\ X'i = v'a^' + a'v\ (6) 
where w% a* are the vector and axial- vector coupling constants for the photon and Z^: 

= 1, a'' = 0, (7) 

= (-l + 4s^)/(4sH^CH/) 

= -l/iAswcw) (8) 

and sw {cw) is sine (cosine) of the Weinberg angle. 

It is convenient to rewrite the cross section in terms of four Born-level matrix elements Mf: 



-Re^(M° + M°), (9) 



where matrix elements are obtained from by crossing symmetry t ^ u. The matrix elements Mj^ are expressed 
through the "even" and "odd" functions and Mo'- 

Ml = £)^*(DT*Mf'^^''-L''^"MfT^T), 

= Dzt^^jjztj^zzzz _ j^zuj^jzzzz^ _ ^^^^ 
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The matrix elements Me and Mg are defined in such a way that they can be used, with minimal modifications, in 
both Born and first-order matrix elements. They are defined in terms of the couplings A 

jkl _ 0('o2 I „,2^|\ ij \ kl.r)/2 ,,2\ -i ij \ kl 



= 2(,S^+7/^)AigAi^'+2(s^-W^)A2'iA2?.', (11) 
ij \ kl . ^ ij -i kl N 



Mvki ^ -2s\X,%X,'^ + X2%X2'^) . (12) 



The kinematic variable y is defined as 

t 1 - cos e i;*' 

y = —,^^^^' (13) 

where & is the center of mass (CMS) scattering angle of the detected electron with momentum k2- E*{E* ) is the 
energy of the initial (detected) electron in CMS, respectively. In the Born approximation, the M0ller scattering is 
elastic and E* = E* , therefore 

1-cose e' , ^, 

yBorn = ^ = 1 - (14) 

where E and E are the initial and scattered electron energies in the Lab frame, respectively. Whenever possible, we 
ignore the electron mass m (which cannot be done in the coUinear singularity regions). 

At the Born level, the unpolarized (averaged over helicity states) cross section is given analytically by 

,^7'., (l + / + (l-y)^) . (15) 

sr(i -y) 

Polarization asymmetry A^n is conventionally defined as 

. <^LL + <^LR — CTRL — CTRR 

Alr = , (16) 

ctll + cflr + (^rl + (trr 

where the first helicity index refers to the beam electrons, and the second helicity index corresponds to the target 
electrons. Since the target helicity is summed over, defined by Eq. p6|l is a parity- violating observable ^2^. The 
Born-level asymmetry is given by Q 



(17) 



where A'^{Q^,y) is an experimental acceptance-dependent analyzing power. Kinematics of the E158 experiment 
correspond to the laboratory beam energies of i? = 45(48) GeV and CMS scatteringangles —0.5 < cos 8 < 0, or 
average invariants s = 2mE « 0.048 GeVVc^ = -t = 0.026 GeVVc^ and y « 0.6 H. 

We perform our calculation in the on-shell (OS) renormalization scheme, defining the weak mixing angle to all 
orders in perturbation theory as cw = m-w/'m-z, where mw and mz are the physical masses of the and 
bosons, respectively. For consistency with the precision electroweak measurements ^J^, we use 

mw = 80.390 GeV (18) 

mz = 91.188 GeV (19) 

which implies 

sin^ 9w = 0.2228 (20) 

in the on-shell scheme. As we will note in Section IVll while the absolute value of the Born asymmetry depends on the 
electron neutral current coupling 1 — 4 sin^ 9w , the relative corrections to the experimentally measured asymmetries 
computed here are quite insensitive to the choice of couplings or the renormalization scheme. 



III. RADIATIVE CORRECTIONS 



It is well known that effects of "internal" bremsstrahlung (real photon emission) need to be combined with the 
contributions from the other one-loop (leading order) electroweak radiative corrections (so-called virtual, or V- 
contributions). The virtual contributions to M0ller scattering have been studied extensively |^ 0] and are not 
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repeated here. However, we need the infrared-divergent part of the V-contributions, which cancels out the corre- 
sponding divergences in the bremsstrahlung cross section. NaturaUy, the extraction of the IR-divergent parts is a 
somewhat subjective procedure, which may lead to ambiguities similar to the concept of scheme dependence in the 
ultraviolet renormalization. We follow the framework of Bardin and Shumeiko jTll | and extract the infrared-divergent 
contributions that are strictly proportional to the Born cross section. Such contributions cancel in the parity- violating 
asymmetry, and remain small even after other corrections are taken into account. We take the infrared-divergent parts 
of the virtual corrections from |1(t1 I17|. which include the vacuum polarization and vertex correction contributions. 
We also compute the IR-divergent 77 and box diagrams. 

The virtual contributions to the M0ller scattering cross section and asymmetry can be classified into three categories: 
the vacuum polarizations of the gauge boson propagators, vertex corrections, and box diagrams (see Fig.Hl): 

= + a^'^"'' + . (21) 



The vacuum polarizations of 7 and Z'^ bosons are infrared-convergent, and arc not considered in this paper. The 
infrared-divergent parts of the vertex corrections (Figures and Et) are obtained from form- factors SFyj^ given in 
[T^ (for fc^ = t, u). Substituting the coupling constants for the vertex form-factors {e.g. v"' SFy'^) in the expressions 
for the Born functions Mg^o, we get the vertex part of the cross section 

a^^^ = ^Re^(A/,^ + 0, (22) 
1=1 

where 



MY = i?^*(Z?^*(Afr^^^ + A4-'^^"'^)-i?^"(Mr^^^ + AfJ^^''^)), 

MY = D^\D^\Mf^^^ + M^^P"^)-D^^{Mf^^^ + M^^P"^)). (23) 

The box diagrams with at least one photon (i.e. Figures |3i andl^t plus u-channel graphs) also contain infrared 
divergences. The diagrams with two Z° or two W bosons are infrared-convergent and are not considered here. We 
compute the box diagram contribution as a sum of the infrared-divergent and infrared- finite parts = cr|? -I- (rfj^. 
The IR-finite parts of the 77 and 7Z boxes are expressed by 

= E {B^n + B'.z) + {t^u). (24) 

k=l,Z 
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The terms B have the following form: 



^(77) 



ilZ) 



kt 



+ (77) 



At low energies (s, |u| ^ J^^^J^ are given by the following expressions: 



°(77) 



°(77) 



L 



2 + 



2t 



2t 



The logs in the electromagnetic box diagrams are 



Ls = -^a; = In^ 



and the scalar integrals in 7Z-parts are 

IjZ = 
IjZ — 

P = 



1 



zdz I dx — = In 



XZy^ — U + vT? 



4(1 - z)(tz{x- 1) + m|) 



(25) 



(26) 



(27) 



(28) 



Over a fairly broad kinematic range of interest lO"'^ GeV^ < \t\ < 1 GeV^ Eq. ^ can be integrated numerically 
and approximated by the following expression to better than 1% precision: 



i-yZ 



0.123 / , 

1.64 + In ^ 



— 3- 0.102 + In ^ 

V —u 



Separating the infrared-divergent and finite virtual contributions, we can write 

a^ = ar« + a^(A^-.), 



(29) 



(30) 



where we introduce the finite photon mass A to regulate the IR divergence. The IR-divergent part proportional to 
the Born cross section is 



^IR 



2a s tu 
In -2 (In - l)cr^ 



(31) 



IV. BREMSSTRAHLUNG CONTRIBUTION 



The complete leading order radiative corrections need to include the inelastic processes that correspond to the real 
photon emission e^e^ —^ e~e~7 (/^-contributions, or real photon bremsstrahlung). The diagrams are shown in Fig.^l 
(plus crossed terms for a total of 16 diagrams with 7 and propagators). Let k be the four- momentum of the 
emitted photon. The differential cross section is given by 




where indices i and j refer to a particular contribution to the cross section (u and t channels with Z'^ or 7 exchange): 



«>J = (1;2;3;4) = {-/t;ju;Zt;Zu). 



(33) 
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FIG. 3: t-channel bremsstrahlung contributions to process e e e e ^. M-channel contributions are obtained by crossing. 



We also use a somewhat unconventional set of kinematic variables 

Zl = 2fcfci 

ti = {p2~Vif 
V — 2kp2 = s + u + t — 4m^ 
Vl = 2kpi = s + u + ti — 4to^ 

Z2 = {ki - P2)'^ = U - V + Zl , 



(34) 



which satisfy the following relations: 



Vl — V — Z — Zl — ti — t . 



(35) 



The integration variable v in Eq. (|32|l describes the inelasticity of the reaction, the deviation from the 2-body 
kinematic constraint s + u + t — Arri^ . The kinematically allowed region for variable v is given by 



V < V 



lim 



2(s + t-4m2) 



1 + - 4m2/s)(l + 4m2 /t) 



s + t 



(36) 



The limit v = w'™ corresponds to the coUinear singularity u = 0. However, for E158, which only detects scattered 
electrons with laboratory energy E' > Ecut ~ 11 GeV, the integration region is further restricted to u™^'^ — 2m(m — 
Ecut) and 



s + t + u" 



The squares of the matrix elements are given by 

M^'j = (M«).. + (M«),. + (M«)., + (M«).„ 

for ij — 13, 31, 11, 33, and 

Af« = (M^)/ + (M«)i + (M^),„ + (Af«), 



(37) 



(38) 



(39) 



for ij = 12, 14, 32, 34, where the traces of the appropriate 7-matrix combinations are multiplied by the density matrices 
and the corresponding propagators 
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where 



V 'zz 

(M^)^^ = Tr[GrP^^(fci)7.A(fc2)]Tr[7^p'^(pi)Gr'^A(p2)]i?^*^i^^'*, 

(M^)^^ = Tr[GrP*^'bi)7.A(p2)]Tr[7^p^^(fci)Gr^A(fc2)]i?'*i?^*\ 

(M«)^^ = Tr[GrP*^(pi)Gr''A(p2)]Tr[7^p^-'-(fci)7.A(fc2)]i?'*i?^*, 

= Tr[GrP^^'(fci)GrA(p2)7pP^^(Pi)7.A(fc2)]i?^*^i^^", 

(M,f), - Tr[GrP^^(fci)7.A(p2)7MP^''bi)GrA(fc2)]i?**^D^'^% 

(M«)^,^ = Tr[7^p^^(fci)G3^"A(p2)Grp'^(pi)7.A(A;2)]i?^*i5^", 

(M«)^ = Tr[7^p^^(fci)7,A(p2)GrP^^'(pi)GrA(fc2)]i?^*i?^"^% (40) 



A(p) ^p + m, p = 7%, (41) 



— 7 1 ^ 



-21 



(42) 



Gr =7'^ 



(43) 



Gg = 7 +7 



(44) 
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2p?-7"/c ,2fc£+fc7° 
-7+7 



(45) 



Expressions for for other values of i and j {e.g. {ij} = 22, 44, 24, 42, etc.) can be obtained from the symmetry 
of expressions in Eq. H4U|I : 



fc2^P2,fcl«Pl 



and 



(M2l)„, = (M«3) 



22 It^U ■ 



(46) 



(47) 



Cases {ij} — 22, 44 are analyzed in a similar manner. The symmetry noted above is less apparent in the interference 
terms (indices zv and vz): 



(M2l).„ + (Mf2)., = (Mi1).. + (A^3l) 



(48) 



and likewise for other parts in Eq. 132|) . 



A. Infrared Divergences in Bremsstrahlung Contributions 



Next we need to address the issue of the infrared divergences in the bremsstrahlung cross section. According to the 
prescription of Bardin and Shumeiko , we find the infrared-divergent parts in the squares of the matrix elements 
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that are proportional to the corresponding Born contributions: 

,2 



(Ml 



IB, 
IB 



where 



Mil 

Mi3 

M31 
M33 
M12 

Mi4 

M32 

M34 



u 



4( 



zf ZZl ' 

s s 



u 



(49) 



ZUl 
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TO 

y2 

2u 

Z\V 

2t 

ZlZ 

2u 



ZlVl ZlV 
ViV 



)M, 



I] 1 



Am' 
, 2s 

+ + 

ZiVi - - 

, 2s 



ZlVi 



vt 



ZiV 

2u 

VlZ 



2t 

ZZl 

2u 

ZVi 

2t 

VVl 

2t 

VVl 



ZV 

4to^ , 
4to^ 



V 

2s 



J^ltjjZtj^Z-yZ-y 
j~^Ztj-yZtj^ZZZZ 
jjit £)7"_/\^7777 
jjit jjZu j^,pZ^Z 
Dl^JjZt j^Z-iZ-i 
j-)ZtjjZuj^jZZZZ 



(50) 



The complete i?-contribution to the cross section is infrared-divergent, but can be separated into an infrared-infinite 
B 

IB 



part cr|^ and the finite part Cp 



'IB 



(51) 



The IR-divergent part of the bremsstrahlung matrix elements, proportional to the Born contributions, can be con- 
structed from Eq. (|49|l according to Eq. I|32|l . The finite contribution to the cross section is then obtained by 
subtraction 



(M«)^ = M^-(M«)'^, 



(52) 



The infrared-divergent part of Eq. H51() , integrated over variables k and v is given in terms of a finite photon mass 
A by 20] 



'IB 



2a , If 
"^"^ V2 



max) 2 



Da" 



(53) 



The integration over the phase space of the bremsstrahlung photon is performed analytically, and integration over 
variable v is done numerically due to complexity of the integral expressions. The photon phase space integral can be 
written as l21il 



I[A] 



ko 



S[iki 



Pi 



kf-m^][A] 



dti 



dz 



(54) 



where Rz is the Gramm determinant (modulo —1), and can be parameterized as a second-order polynomial in z as 

Rz = AzZ^ + 2BzZ + Cz. (55) 
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Coefficients A^, B^, and Cz are given by the following expressions: 

Az = -im^t- {v~tf, 



Bz = Eti+F 

E v(u - 2m^) ~ st, F ^ t{2m'^v + st ~ sv), 

C, = -{Atl + 2Bh+C), (56) 

A = {s- v)^ - 4m^s , 

B = —st{s — w — 4m^) — 2m^v^ , 

C = si2(s-4m2) . 



Integration limits z™'" ^ and i™'" ^ '"'^'^ are solutions of the equations Rz = and z" 



^min/max ^ ± \/ B^ - AzCz ^^^^ 

Az 

min/max _ - t^) + 2mH =F f V-^z 

2[v + m/) 

This set of variables makes integration more convenient. 

Expressions for (M^) are computed using symbolic manipulation program REDU CE j22l| but are excessively complex 
to be listed here. They are available as subroutines in the Fortran program rcAPV[2^. The relevant integrals, listed 
in the Appendix of Ref. , were computed both analytically and numerically. 



V. CANCELLATION OF INFRARED SINGULARITIES 



Adding infrared-divergent parts of V- and i?-contributions given in Eq. and Eq. H53|ll together with the IR- finite 
pieces, we obtain the finite expression for the radiatively-corrected cross section, free of non-physical parameters: 



where 



<^IR + ^IR + ^F+^F + (59) 

^ (4 In ^ (In ^ ~ 1) -f 5f -I- )a° + aj^^ + + a|? , 



5f 



£(^_1 s\s + tf 1 -f-^-l + l 



m** 2 
2 



dv 



-tm^ 2 
2 , , 2 



( , , V ^ 'I , , V ^ 'A , , \ , , V ^ 
In 1 - - + - In 1 - - + - In 1 r - - In 1 + ^ 



H ■ In In ^ ^ In ^ ' 



s + t — V 
l„, = In — - . 



s — V m'^T V — t m'^T 



.s + t 2 

In , T = V + m , 



(60) 



VI. RESULTS AND DISCUSSION 



A. Numerical Results 



In the following, we evaluate the effect of the bremsstrahlung radiative corrections on the parity- violating asymmetry 
A]^j^ in the scattering of the longitudinally polarized electrons off unpolarized target electrons. We consider the 
kinematic conditions that correspond to the experimental setup of the SLAC E158 experiment p^ . i.e. beam energies 
of 45 and 48 GeV. E158 setup is such that the radiated photon is not detected. Moreover, the scattered electrons are 
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FIG. 4: Differential cross section (left) and asymmetry Alr (right) as a function of y for beam energy E = 45 GeV. Dotted 
line shows the Born cross section and asymmetry, effects of bremsstrahlung corrections are displayed by dashed lines, and the 
solid line is the result of the full QED corrections. The plots are cut off at y « 0.75 by the E158 spectrometer acceptance 
(_B' > 11 GeV). 
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FIG. 5: M0ller scattering cross section (left) and asymmetry Alr (right) as a function of scattered electron energy E for 
beam energy E — 45 GeV and y — 0.5. 



only detected if their energy is above the threshold £" > 11 GeV; this restriction limits the range of integration over 
variable v as given in Eq. H37|) . 

The relative corrections to the cross section and asymmetry can be defined as 



6a 



6A 



A 



LR 



40 



(61) 



where Aq is the Born asymmetry, and A^'~' is the radiatively corrected asymmetry. 

The Born cross section a''\ the radiatively corrected cross sections 0-^''°"^ (which includes only soft and hard 
bremsstrahlung corrections and results of IR cancellation) and a^^^ (full QED corrections), as well as the asymmetry 
Alh are shown as a function of variable y in Fig.^jfor beam energy oi E — GeV. Fig.|Slshows the double-differential 
cross section (Pa/dE' /dy and the asymmetry A^n as a function of the scattered electron energy in the lab frame E' 
for a fixed y = 0.5. This double-differential cross section is used to properly average the radiative corrections over 
the experimental acceptance. The corrections to cross section and asymmetry are shown in Fig. The numerical 
precision of the corrections is about 0.1%. 

We find that at fixed value of ?/ = 0.5 and E — 45 GeV {Q^ = 0.023 GeV^), hard bremsstrahlung reduces the value 
of the parity-violating asymmetry by « 1%. Contribution from the 77 and 'yZ box diagrams is also negative and 
reduces the asymmetry by « 5%, so that the total QED correction at j/ = 0.5 and E = 45 GeV is —6.2%. 

At fixed scattering angles, radiative effects move events towards lower values of y. Therefore, even though at a 
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FIG. 6: Relative corrections to the M0ller scattering cross section (left) and to the parity-violating asymmetry Alr (right) 
as a function of y. Dashed line includes bremsstrahlung corrections, and the solid line is the result of the full first-order QED 
calculation. 



fixed value of y — 0.5 the hard bremsstrahlung correction is negative, the net change of the asymmetry, integrated 
over E158 acceptance -0.5 < cos 6 < at S = 45 GeV is {SA^'-'"'^) « +1%. The full QED correction, including 77 
and boxes, is (M^^^) = -4.5%. 



B. Factorization of QED Radiative Corrections and NLO Uncertainties 

At leading order, corrections to the parity- violating asymmetry from the diagrams involving photons (i.e. soft and 
hard bremsstrahlung, 77 and jZ boxes, and photonic vertex diagrams) are proportional to the product v^a^ . In the 
OS scheme, these contributions are strictly proportional to the Born asymmetry. In other words, in the OS scheme, 
the relative corrections to the asymmetry due to QED diagrams are independent of the weak mixing angle, and can 
be factorized out. We can write 

A^l^ = A'{Q',y){l + SA^^^) (l -Asm' Ow) (62) 

where A'^{Q'^, y) is the Born analyzing power defined in Eq. (|17|) and SA*^^^ is the relative radiative correction. For 
E = A5 GeV and y = 0.5, full QED correction is SA'^^^ — —0.061, and the average over E158 kinematics (beam 
energies of 45 and 48 GeV and -0.5 < cos© < 0) is {SAQ^°)ei58 = -0.043. 

Virtual corrections, such as vacuum polarization and box and vertex diagrams with heavy gauge bosons, are not 
in general proportional to v^a^ , and spoil the simple factorization of Eq. (|62|l . Nevertheless, one can relate the 
leading-order asymmetry to the Born-level formula through the set of multiplicative and additive corrections: 



ALR{Q^,y) = ^°(Q^ zy)p(Q2)(l + SA{Q',y)) (l - 4sin2 0^(Q2) + a{Q')) . (63) 

Here the effective mixing angle relevant for a scattering process at momentum transfer Q is defined through the 
form-factor k{Q') (in OS scheme) or k{Q' , fi^) (in MS scheme): 

sin2 0^(Q2) = K(Q2)gj^2^os ^ k{Q'^fj.')smHwi^i') ■ (64) 

In MS scheme, one typically chooses /i = mz- 

All corrections in Eq. (|63|l are of order 0{a) but have a different physical meaning. The form- factor k = 1 + 0{a) 
defines the momentum dependence (running) of the effective weak mixing angle. We would like to define it in a process- 
independent way, such that various experimental measurements could be directly compared in terms of sin^0^. 
Typically, k includes contributions from j — Z mixing and anapole moment diagrams (Fig. |2K-c) , but may include 
other terms [Ullil. In definition of Ref. carried out in MS scheme with sin^ (to|) = s'^ = 0.23120±0.00015 0, 
k(Q2 = 0, = ml) = 1.0298 ± 0.0026 

The form-factor p{Q') = 1 + 0{a) is a low-energy ratio of the neutral weak coupling to the charge coupling. It 
depends on the choice of the Fermi constant Gp in Eq. H17|l : for Gp derived from the muon-decay constant G^, 
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this correction is |^ p = 1.0012 ± 0.0005. p contains logarithmic dependence on the Higgs mass (we use mn = 
1 13^40 GeV 0) and linear dependence on the top quark mass (mt = 177 ±4 GeV 0). 

The remaining first-order corrections are included in Eq. Ht)3|) as factors 5A{Q'^,y) and A((5^). A((3^) typically 
includes box diagrams with two heavy bosons, and 5A{Q^,y) contains the kinematics-dependent factorizable QED 
corrections computed in Section IVII 

The only remaining question is evaluating the next-to-leading order correction uncertainties. Normally, NLO 
corrections would be of order 0{a) of the LO terms. However, we have to pay special attention to the logarithmically- 
enhanced contributions in the 7Z box diagrams, e.g. terms proportional to ln(m2/s) in Eq. (|29|1 . 

The effective Z°-electron coupling, Qw = — p(l — 4sin^ 9w) changes by about 40% between zero momentum transfer 
and Z° pole 0. Since the box diagrams in Fig. |21 involve integration over internal momenta of the photon and 
propagators, a complete calculation of the box diagrams has to take into account the momentum dependence of Qw 
and a. Strictly speaking, the momentum dependence of the weak or electromagnetic charges in the box integrals is 
a next-to-leading order effect, and is beyond the scope of this work. However, a judicious choice of the coupling at 
leading order could reduce the NLO corrections. 

Czarnecki and Marciano @ have argued that the NLO errors are reduced if the 7Z box diagrams are evaluated in 
MS scheme with the average value of the weak mixing angle 

(sin^ 9^) ^ ^■^^'^'Ow{Q)+^mHMml) ^^^^ 

since the relevant integrals over the internal momentum k are dominated by the poles near ^ and k^ ^ m^. 
Similar arguments apply to the value of the fine structure constant in the 7Z box diagrams. Thus, following the spirit 
of this argument, we move the leading logarithmic contribution to the jZ box diagrams j5j] 

A,z(Q^) - ^ (1 - 4sin2 9wiQ')) In ^ (66) 
3 An s 

from the multiplicative correction SA{Q^,y) to the additive correction A{Q'^). Moreover, we use the average value 
{A.yz) = {A.yz{0) + Ajz{'mz))/2 in the expression for A{Q^), and treat the spread \A^z{0) — A~fz{'mz)\/2 as 
an estimate of the NLO uncertainties. This choice keeps all of the experimental acceptance dependence in the 
multiplicative correction 1 -I- SA(Q^ ,y), properly propagates the bremsstrahlung corrections, and reduces the overall 
size of the theoretical uncertainties. For experimental kinematics of the E158 experiment, we find 

M(Q2^2/) = 0.006 ±0.005 , (67) 

where the uncertainty is dominated by the possible variations of the experimental acceptance and numerical precision. 
Employing the calculation of Ref. for the terms k and p, the corresponding contributions from the WW and ZZ 
box diagrams, and our definition of (A^z) term, the residual additive correction in Eq. H63|) is 

A(Q2) ^ (A,z) + - ^"i?^^, (1 - 4g^) [1 + (1 - 4^')] - -0-0007 ± 0-0009 (68) 

47rs^ 327rs^(l — s^) '■ ■' 

where the uncertainty is dominated by our conservative estimate of the NLO terms. 

The fact that both corrections SA{Q'^,y) and A{Q'^) are small is somewhat accidental. SA{Q^,y) is dominated by 
the bremsstrahlung contributions and the infrared parts of the 77 and 7Z diagrams that happen to cancel each other 
for the asymmetric acceptance of the E158 spectrometer. On the other hand, A{Q'^) is small at the of E158 due to 
the cancellation between the WW box diagrams and the large logarithmic contribution in the "/Z box. Such precise 
cancellation is not expected at lower momentum transfers, or for other processes, such as elastic ep scattering. For 
example, at y = 0.5 and E = 12 GeV, which corresponds to the idealized kinematics for a proposed M0ller scattering 
experiments at the Jefferson Lab |2^, the corrections would be 

M(g^ J/ = 0.5)12 Gov = -0.024 ± 0.005 , A{Q^)i2 GcV = -0.0011 ± 0.0010 . (69) 

For a 250 GeV fixed target M0ller experiment, e.g. at a future Linear Collider |23|, one would find 

5A{Q^,y = 0.5)250 GcV = -0.012 ± 0.005 , A{Q^)i2 GoV = -0.0002 ± 0.0007 . (70) 



VII. CONCLUSIONS 



In conclusion, we have computed the QED corrections to the parity- violating left-right asymmetry A^n in M0ller 
scattering. We used a covariant method for removing infrared divergences without introducing unphysical cutoffs. 
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For the kinematics of SLAC E158 experiment, the overall corrections appear to be small, due to a fortuitous cancel- 
lation between electroweak and electromagnetic terms. We reduce the theoretical uncertainties due to higher order 
logarithmic terms by the appropriate choice of the couplings used to compute the box diagram contributions. Our 
calculation is applicable to a wide range of fixed target energies, from the proposed M0ller scattering experiments at 
12 GeV at the Jefferson Lab (26(1, to the possible fixed target experiments at a future Linear Collider |27l |. 
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